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Abstract 

. In a previous paper we studied the collapse of a spherically symmetric dust distribution (marginally 

bound LTB) in d-dimensional AdS spacetime and obtained the condition for the formation of trapped 
surfaces. Here we extend the analysis by giving the canonical theory for the same and subsequently 
quantize the system by solving the Whcclcr-DeWitt equation. We show that for the case of small dust 
perturbations around a black hole the wave functionals so obtained describe an AdS-Schwarzschild black 
' hole in equilibrium with a thermal bath at Hawking temperature and show the non-trivial dependence 

of this temperature on the number of spacetime dimensions and the cosmological constant. 
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1 Introduction 



One of the fundamental unresolved questions in physics is whether there exists a consistent formulation of 
quantum theory of gravity and, if it does, whether it cures the singularity problem that plagues the classical 
theory. Also it is hoped that such a theory, if it exists, will give fundamental insights into the nature of 
Hawking radiation pQ and the origin of black hole entropy as predicted by the use of semiclassical methods 
■ 0) [2]) @]- As long as a fully consistent theory of quantum gravity is not available it is worthwhile to 
study simplified models using different techniques. Among the most useful models are those with spherical 
symmetry (midisuperspace) , since the assumption of spherical symmetry keeps the analysis at a technically 
simpler level and yet allows one to study various non-trivial issues that are expected to be important in 
O quantum gravity. 



The midisuperspace quantization program has been carried out in a series of papers for spherical dust 
collapse (LTB collapse) in 3+1 dimensions, both for the marginal [5], [6], [7], [8], [9] and the non-marginal 
case |10| . following the work of Kuchar for the Schwarzschild black hole [12]. There one could set up 
and solve the Wheeler-DeWitt equation for these models and thus obtain exact quantum states, which in 
the near horizon approximation gave Hawking radiation. However, issues relating to singularity avoidance 
and to the exact quantum gravitational nature of Hawking radiation could not be addressed fully due to 
various technical difficulties. To overcome some of the technical difficulties that arise in 3+1 dimensions, 
the quantum gravitational dust collapse model has also been studied in 2+1 dimensions [13] where it 
turns out that a black hole solution (the BTZ black hole) is only possible in the presence of a negative 
cosmological constant. However, in this case the thermodynamics and the statistics of the quantized BTZ 
black hole are completely different from that of the 3+1 dimensional Schwarzschild black hole. Specifically, 
whereas the Schwarzschild black hole follows the Boltzmann statistics, the BTZ black hole follows Bose- 
Einstein statistics [14]. This is possibly related to the fact that the BTZ black hole has positive specific 
heat, whereas the specific heat is negative for the Schwarzschild black hole. 

Motivated by this observation we were led to the question: what is the role of the cosmological constant 
and the number of spatial dimensions in determining the nature of thermodynamics and the statistics of the 
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quantized black hole? The present paper is the second in a series of three papers addressing this question. 
In the first paper [15] we solved the Einstein equations for a collapsing dust ball in an asymptotically 
AdS spacetime and examined the nature of the gravitational singularity. Here we look at the canonical 
formulation of the same which we then quantize obtaining the Wheeler-DeWitt equation, which is then 
solved on a lattice to obtain exact quantum states. These quantum states are shown to give rise to Hawking 
radiation in the limit of small dust perturbation around an AdS-Schwarzschild black hole. We also allude 
to the interplay between the number of spatial dimensions and the value of the cosmological constant and 
the black hole mass in determining the thermodynamical properties of the black hole. The statistical origin 
of this dependence based on the present quantum treatment will be taken up in a subsequent paper. 

The paper is organized as follows. In Sec. II we recall the important results concerning the collapse of 
an inhomogeneous dust cloud under the assumption of spherical symmetry in d = n + 2 dimensions in the 
presence of a negative cosmological constant. The canonical formulation for this system is discussed in Sec. 
Ill where we also give the contribution(s) to the action coming from the boundary variations along with a 
discussion of the fall-off conditions. In Sec. IV we make a canonical transformation following Kuchar to 
obtain a much simpler expression for the Hamiltonian constraint, which is followed by a discussion of the 
dust action in Sec. V. In Sec. VI we present the quantization of this system which leads to the Wheeler- 
DeWitt equation. This equation is then solved to obtain exact quantum states on a lattice following 
DeWitt's regularization scheme. Subsequently we obtain Hawking radiation from these states using a 
particular ansatz for the inner product. It is remarkable that the d-dimensional spherically symmetric 
system with a negative cosmological constant can be cast in the same canonical form as the 3+1 system 
without a cosmological constant. 



2 The Classical Solution 

The metric for a spherically symmetric inhomogeneous dust cloud, described by the stress-energy tensor 
e(r, p)UuJJ v , in d-dimensional spacetime (d = n + 2), in the presence of a negative cosmological 



constant —A (A > 0), in co-moving coordinates is 



ds 2 = 



-dr 2 + 



R 2 



1 + 2E 



dp 2 + R 2 dn 2 n 



(1) 



Here f2 n is the solid angle for the n-sphere, R = R(t, p) is the radius of the n-sphere and E, the so- 
called energy function, is some function of the radial coordinate p and we use the notation R = d p R and 
R* = d T R. For the marginally bound case, corresponding to E = 0, Einstein equations imply that the 
energy density e of the dust cloud and the metric component R evolve as (see |15j for details of the classical 
solution) 



P) 
(R*) 2 



(n - 1) F 
8vrG RnR 
2A o2 



n 



(n + V 



R 2 + 



F(r) 



R< 



(2) 
(3) 



Here F, the mass function, is some function of p with the condition F > 0. We note that collapse is 
described by R* < 0. Using the scaling freedom in the choice of the radial coordinate p to set R = p at 
the initial time r = we can write the solution of the last equation as 
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Equation @ shows that for a shell labelled p the curvature radius R becomes zero at the time 



sm -1 . / JAp_"+i 



T o(p) = 7== • (5) 

' A(n+1) 



2n 

This solution can be matched to an exterior Schwarzschild-AdS spacetime 

, ■? ( F(p b ) 2Ax 2 \ , 9 / F( Pb ) 2Ax 2 \ _1 , 9 9 _ 9 

ds 2 = - 1 ^ + — WiT 2 + 1 ^ + — dx 2 + x 2 dfl 2 (6 

V a;"- 1 n(n + 1)/ \ x"- 1 n(n + 1) J w 

where (T, x) are the coordinates in the exterior and p^ is the boundary of the dust cloud. 

3 Hamiltonian Formalism 

We now want to set up the canonical theory for the general d = n + 2 dimensional, spherically symmetric 
collapse problem with and without the cosmological constant. For this we consider the ADM metric 

ds 2 = -N 2 dt 2 + L 2 (dr + N r dt) 2 + R 2 dn 2 n (7) 

where as before £l n is the solid angle 

dO, 2 n = ddf + sin 2 6i{d6 2 2 + sin 2 9 2 (d9 2 . . . sin 2 0„_i(i# 2 ) . . .)). (8) 

The Einstein-Hilbert action in presence of a negative cosmological constant —A (A > 0), (ignoring the 
boundary terms), is 

Ser = [ d d x^( d K + 2A) (9) 



where d lZ is the Ricci Scalar of the d-dimensional manifold and Gd is the gravitation constant in d- 
dimensions. The above equation after using and integrating over the angular dimensions becomes 



n-l 



Sek = ^ r( ^ } ^ J dt J drNLR n ( d - l U + K^K^ - K 2 + 2A) = J dt J dr C (10) 

where d ~ l lZ is the Ricci Scalar of the spatial slice and the prefactor comes from integrating over the solid 
angle. K^ u is the extrinsic curvature. If we define 

r \ ri—l n+1 

~! _ (n - 1)tt 2 _ 2vr 2 

U d - 4nr( n±l )Gd < iZ «- r (n±l) 



then 

C = 

Gd 



C = i NLR n { d - 1 TZ + K^ V K^ V -K 2 + 2A) (12) 



or 



nR n ~ 2 

C = — ^ (A^T?' - J R)[2(AT r Li?) , + (n-3)7V r Li? / -2i?L- (n- l)Li?] 

at nn—2 

^ [-n(n - 1)L 3 - 2nRL'R' + n(n - l)LR' 2 + 2nLRR" - 2KR 2 L i ] (13) 

GdL 2 
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From here, we compute the momentum densities, 

2nR n ~ 2 



Pr 
Pl 



G d N 
2nR n - 1 



G d N 

which can be solved for the velocities, 



[(N r LR)' + (n - 2)N r LR' - RL - (n - 1)LR], 
(N r R' - R), 



(14) 



R = -^4^+N^R', 



2nR 
G d N 



Pr ,, U LP L 



+ (N r L)'. 



2nR n ~ l 2nR n 

With these definitions the Lagrangian density can be expressed in terms of the canonical variables 



(15) 
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[-n(n - l)L s - 2nRL'R' + n(n - l)LR' 2 + 2nLRR" - 2AR 2 L S ] 



and putting it all together, we compute the Hamiltonian 

H = J dr[RP R + LP i - C] = J dr[NH + N r H r ], 
from which we can read off the Hamiltonain and the diffeomorphism constraints, respectively, 
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It can be shown that this gives the same results for the n = 2 case as we have derived in our earlier paper. 
3.1 Boundary Variations 

Our next task is to list all possible contributions to the boundary variations. From the diffeomorphism 
constraint we obtain, 

- / dtN r P R 5R 

+ / dtN r L5P L (19) 

JdT, 

and, from the Hamiltonian constraint we get, 

+ 2n / dt—^ — -L'SR 



+ 2n / dt 
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The last boundary variation can be written as the sum of three parts: 



2n / dt 
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(21) 



The middle term cancels the third boundary variation of the second list and the last term cancels the first 
boundary variation (in the same list). We are then left with just four potential contributions from the 
boundary, viz. 



2nNR "^ gSL + N*LSP L + (-N*P R + 2 " J !' R "" I 8R 



+ I dt 

Ian I G d L 2 V G d L 

Whether or not any term contributes depends on the asymptotic conditions we will shortly impose. 



(22) 



3.2 Equations of Motion 

We now use the equations of motion q = {q, M}, where q = (R, L, Pr, Pl), to obtain the evolution equations. 
Making the choice N = 1 and TV 1 " = we obtain 
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Taking the time derivative of (J23]) and substituting for Pl we find 

R- 



R 2 n-1 (n - l)R' 2 AR 
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(23) 
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We need to do a similar thing for L and so we consider 



?L = GdP L +( n - l) GdPh 
R' 2nR n ~ 1 R' { } 2nR n ' 



Using the momentum constraint Tt r ~ to eliminate P' L we obtain 

B! G d P R G d P L L 

T7 + ( n - 1 r, 



R' 2nR n ~ 1 L 



2nR n L 



(28) 



(29) 



This implies L = y- where / is a function only of r. If we take f(r) = \Jl + 2E(r) and substitute the 
expression for L in (|27[) we find 



- . . R? , . E AR 
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v 1 2R y ' R n 



This equation can be integrated with respect to time 
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This is seen to be the correct equation for the evolution of R in the non-marginal case. F(r) is of course 
the mass function and E{r) is the energy function (which is zero for the marginal case). 

3.3 Fall-off conditions 

To determine the fall-off conditions at infinity, we imagine that the metric is smoothly matched to the 
Schwarzschild-AdS metric at some boundary given by 



where 



ds 2 = -f(r)dt 2 + f- l (r)dr 2 + r 2 dn 2 n 
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where we have set F(rb) = 2M = G^Mjn and where M is the ADM mass. The fall-off conditions at 
infinity are given by (n > 1) 
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Let's now look at the boundary terms in (|22p . Inserting the dependences, we see that the contributions 
amount to 

^ f dtN + (5M+ , 2A 6^) . (35) 
G d Jan V n(n + 1) J 



We will call 



2n 

M + = 



77 2AR l 
M + 



(36) 



nyn + 1) 

so that the contribution at infinity is 

dtN + 5M + (37) 



IdT, 

and must be subtracted from the hypersurface action. 

The fall off conditions at r = depend sensitively on the choices we make for the arbitrary functions 
E(r) and F{r). Near r = 0, take E(r) = Yl^Lo ^"nX n and F(r) = Yl^Lo ^i" n , requiring Fq > (avoiding 
shell-crossing singularities etc. will involve appropriate conditions on the coefficients). In this way we 
determine the following fall-off conditions to be met: 

R{t,r) -^a{t) + b{t)r + O (r 2 ) 
L(t,r)^ 1 b(t) + O (r) 
P R (t,r)^P R0 (t)+O°(r) 
P L (t,r)^P L0 (t) + O°(r) 

N r (t,r) -» O (r) (38) 

where 7 = — 2Eq and is thus equal to one for the marginally bound case we are considering. They 

guarantee that the hypersurface action is well behaved at r = 0. Inserting these into (j22[) gives 

-Jr- / dtN 5lnb. (39) 

If we call 5 In 6 = —G d 5Mo/2n then we must add the term 

dtN 5M (40) 

as 

to the hypersurface action to cancel the contribution from the origin. Combining the two, we must add 
the term 

Sbe = - [ dtN + M + + f dtN M (41) 
to the hypersurface action to cancel unwanted boundary variations. 

4 Canonical Transformations 

Our aim now is to write the mass function F in terms of the canonical variables. For this we start by 
embedding the ADM metric (|7|) in the spacetime described by the LTB metric JT]). In what follows we will 
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use the prime and an over-dot to refer to derivatives with respect to the ADM labels r and t respectively 
and we define R = R/y/1 + 2E. With this a foliation described by r(t,r) and p(t,r) leads to 



L 

N r 
N 

Using these expressions in the momentum densities in (|14p we find 
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R*i + Rp = R*f + Ry/1 + 2E p, 
R> = R* T ' + Rp' = R*t' + Ry/1 + 2E p' . 
Using this along with (|31 1) in (|45|) we find that 
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where the upper sign is appropriate for expansion and the lower for collapse. Solving for r' and calling 

F 2AR 2 , t . 
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=F-RVl + 2E-F- 
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This is a generalization of the result for n = 2 in [10] to arbitrary n. 

If we now substitute this in (I42p . taking care to maintain the appropriate signs for expansion and 
collapse, we find for either case 



Solving for T we have 
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and therefore the mass function is obtained in terms of the canonical variables as 
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We would now like to turn the mass function into a canonical variable. Taking the appropriate Poisson 
brackets shows that 

* = w (53) 
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is canonically conjugate to F (i.e. {F, P f }pb = 1)> which can be seen to be a generalization of the result 
in earlier work. 

We are looking for a coordinate transformation that would take our diffeomorphism constraint to 



H r = R'P R - LP' L -» R'P R + F'P F 
i.e., we require a new momentum Pr such that 

LP' L F'P f 
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and the set {R, Pr, F, P f } forms a canonical chart on the phase space. We can simplify the expression 
above to get 
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We may check this expression with the expressions we had in 2+1 and 3+1 dimensions. 
• For n = 1 we find 
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For n = 2 we have 
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Both these results were obtained earlier |13j . jlQ] . 

We now set about showing that the transformation is canonical. Since we already know two coordinates 
and one conjugate momentum, we simply use the analogues of the standard equations: 
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(where the left hand side are the old phase space variables) i.e., the equations 
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(63) 



and ask for a generating function ¥[L,R,Pl,Pr]. The last equation in (j63l) tells us that F = W[R,L,Pl\. 
The third equation gives 
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Using this in the first equation (|63[) then implies 



~5L 







(65) 



(66) 



showing that Fi = Fi[i2]. This is a trivial dependence of since we are interested in computing Pr. We 
thus take Fi to be independent of R, a constant. Calculating Pr in ([63]) using the form for F in ([65]) we 
find that we get the same result as we had before in (|60|) . Thus we have established that {R, F, Pr, Pp} 
form a canonical chart. 

Finally we want to rewrite the Hamiltonian constraint in (1181) in terms of the new variables. Eliminating 
L and Pl in favour of F and Pf we get 
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where it is understood that L is expressed in terms of the new canonical variables. We note that it is 
remarkable that the constraints have the same form as the constraints in 2+1 dimensions and in 3+1 
dimensions. 



4.1 Boundary Action 

Because varying N+ and iVo in (I4ip would lead to zero ADM mass and restrict Fq to zero, respectively, 
both iV + and iVo should be considered as prescribed functions. By the fall-off conditions, the lapse function, 
-/V r , is required to vanish both at the center as well as at infinity. This implies that the time evolution is 
generated along the world lines of observers with r = constant. If we introduce the proper time of these 
observers as a new variable, we can express the lapse function in the form N + {t) = r + and Nq{€) = tq. 
This leads to 

S ds = - J dtM+{t)T+ + J dtM (t)t . (68) 

Thus we remove the need to fix the lapse function at the boundaries. Extending the treatment in [12J, 
the aim is to cast the homogeneous part of the action into Liouville form, finding a transformation to new 
variables that absorb the boundary terms. This can be done by introducing the mass density T = nF' /Gd 
as a new canonical variable. Define 

F(r) = —M + — f dr'Tir), (69) 
n n J 
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and reconsider the Liouville form 



6 := -M + 8t + + M 5t + / drP F 5F 
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we can write the Liouville form as 
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where 



Po = -tq, 

p + = t+ H / drP F (r), 

n Jo 



4 See Kuchaf [12]. Consider 



dr'Sr(r) x / dr'P F {r) 



SF(r)x J dr'P F (r') - P F (r) x y dr'<5r(r') 
Integrating the left hand side from to oo gives zero, therefore 
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P r (r) = -— f dr'P F {r'). (75) 
n Jo 



The new form of the action is then 



Ser = J dt (p M + p + M + + J dr [P R R + P T t - NH 9 - N r H r ^j , (76) 
where the new constraints read 

W = -- [-tp r p v + f- x b!t] ps 0, 

H r = R'P R - TPf ps 0. (77) 

Notice that Gd (and n) dependence has been absorbed in the definitions of T and Pp. In this way we arrive 
at the remarkable result that the constraints have the same form in any dimension. So far we have not 
included the dust action. We do this below. 

5 Dust Action 

In this section, we rewrite the Hamiltonian constraint in quadratic form including the dust action. For 
this we begin with the dust action (see [16], [17] for details on the dust action) 

Sd = -\J d d x^~g e(x)[g a pU a UP + 1] = J dt j drNLR n e{x)[g^U a U p + 1]. (78) 
We only consider non-rotating dust and use 

U» = -r )At (79) 
where r is the dust proper time. Using the ADM metric for spherical symmetry we find 



S d = -^- j dt j drNLR n e{t,r) 



r 2 2N r . (I N r2 , l2 



(80) 



and it follows that 

P r = £ = ™^_ NV] . (81) 



Thus 



and the Hamiltonian can be written as 



NP 

+ NW (82) 



n n LR n s 



U = NH d + N r H d (83) 



with 



n d 



H d = t'P t . (84) 

Thus we have the full gravitation+dust constraints 

1 



W - L L 



FP R P[< + F- l R'Y - P t Vt' 2 + L 2 psO, (85) 
H r = t'P t + R'P R -TP^O. (86) 
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We can use the diffeomorphism constraint to simplify the Hamiltonian constraint. Doing so we finally obtain 
the following expression for the Hamiltonian constraint which, because of the absence of any derivative 
terms, is much easier to use for quantization 



n 9 = P 2 + TP\ 



0. 



(87) 



This is precisely the same form that we had in the 3+1 case; thus quantization can proceed as before. 

6 Quantum States and Hawking Radiation 



We now quantize the diffeomorphism constraint (|86j) and the Hamiltonian constraint (187|) . For this we 



make the substitutions P T - 
the diffeomorphism constraint is 



Pr 



-ih 



STjrj 



. With this the quantized form of 



ih 



R'- 



+ r 



6 



^[r(r'),R(r'),T(r')] = 0. 



5R(r) <5r(r) \5T(r 
The quantized version of the Hamiltonian constraint, the Wheeler-DeWitt equation for the system, is 



(88) 



-ih 



5r(r) 



-ih- 



SR(r) 



h 2 A(R, F)6(0) J ^- ) - h 2 B(R, F)5(0) 2 - ^ 



y[T(r'),R(r'),T{r')] = 
(89) 



where the second and the third terms on the left take into account the factor ordering ambiguities in ([8 
We assume that the wave functional is of the form 



y[r(r),R(r),r(r)] = U[ drT{r 



iW(r(r),R(r),T(r)) 



(90) 



Since we would want to put the system on a lattice such that the wave functional can be written as the 
product of wave functions for each lattice point, that is ^ = Y\ ri ipi where i designates the lattice points, 
we take U to be the exponential function. With this choice, and denoting the lattice spacing as a, the 
assumed form of the wave functional becomes 



*[T(r),R(r),r(r)] 



fdrT(r) 



iW(T(r),R(r),r(r)) 
2 



it»^E/ri^ i (T(fi),%),r(r i )) 



lim TT e^^M^)^)^)) 

3 

limn^Kr,),^),^)). 



(91) 



Here we recall that T = nF' /G which implies F(r) is given by (|69p and therefore on the lattice we have 
-F(rj) = lim CT ^o S}=o GaTj/n. We also have 5(0) — > lim CT ^o V "- Now the lattice version of the WDW 
equation is [5] 

' d 2 _ d 2 . , „ „ , d 



V I ^+^ + A(R J ,F,)— + B(R J ,F J )) + 
After substituting for ip~ from (|9ip the above equation gives 



a 2 ^ 



= 0. 



(92) 



+ h 2 Tj 



9Wj 
dR 3 



+ - 



~drf + h:kl ~dRj + h 1 



dRi 



+h 2 B(R j ,F j ) = 
(93) 
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For the above equation to be satisfied independent of the choice of a the following three equations must 
be satisfied f; 10 1 



dW 



Or 

( d 2 



dW 

l ~dR 



4 



0. 



B = 0, 



(94) 

(95) 
(96) 



where the label j has been ignored in writing the above equations, it being understood that these equations 
are to be satisfied at each lattice point. To solve the above equations we look for solutions of the form 
W(r,R) = a{r) + P(R). With this ([94]) and (JHSJ) respectively become 



\dr ) 

d 2 a 



0. 



h T — — 

dr 2 dR 2 



A— = 0. 

dR 



(97) 



(98) 



In both the above equations the first term is dependent only on r while the rest of the terms are independent 
of it (depending only on T and R), which therefore implies that these must be constant. With this (|97p 
implies ti 2 f^) = c 2 where c\ is a constant. The solution of this is a = ±(c\T + cq) / h. Using this solution 
for a, ([97]) becomes 



From this we get 



c\ + n 2 T 



d£ 
dR 



(dp 
\dR 



4 



0. 



This can be differentiated with respect to R to obtain 



d 2 p 
dR 2 



nr 2 



2^4 - c{T 



A-c 2 T 



dT 
dR' 



(99) 
(100) 

(101) 



Also from the solution for a we see that the first term of (j98j) is zero which therefore implies 

d 2 p dp 

Substituting for dP/dR and d 2 P/dR 2 in the above equation we get 

2A(A - c\T) = [c 2 F+ 2(4 - 



(102) 



(103) 



This equation when solved for A gives 



.4 



(n-l)F 



AAR \ 
n(n+l) J 



2c\kR 2 \ 
n(n+l) / 



9^2 , 2c?F _ i4AR?\ 
ZC 1 n(n+l) J 



(104) 



Thus for this A we get a separable solution for W 



W 



cq db c\T ± / dR 



(105) 
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With this form for W the wave function becomes 



/ dr^r(r) [co±ciT±fdR - 

V[T(r),R(r),T(r)]=e V A (106) 

Before proceeding further we would like to express the dust proper time r in terms of the Killing time T. 
From (I49p we recall that (putting E = 0, as is the case for marginally bound model) 

T > = ?L^T- GdLP \ . (107) 

However, we also have Pp = LPl/2R h F = —nP^/Gd, which implies that r can be written as 

r = P rT y dR^y^. (108) 

Now in [T2] it was shown that for the Schwarzschild geometry P? equals the Killing time T (there it is 
2Pp = T; however, recall that we have absorbed n in the definition of Pp). Therefore if we have small dust 
perturbations not affecting the spacetime geometry, this relation may still be used and we get 

T = T^JdR^y^. (109) 

We note that here the minus sign is for expanding dust cloud and the plus sign is for the collapsing dust 
cloud. Now we can also fix the value of the constant c\ appearing in the expression for the wave function. 
For this we note that the Hamiltonian constraint (|87jl can be solved for Pr giving 



P V 2 

PR=± f\h-^ (110) 



This when substitued in the diffeomorphism constraint (|86p gives 



. R' T 2 TPl , , 

Comparing this with (|107p after expressing it in terms of Pp we have that P T = V. This implies that 
P T $? = r^. Evaluating the two sides we get c\ = 2. 
Using this in (|106p the wave function becomes 

fKrj^M.rM^^^^t^N^)^^^]. (112) 

We are interested in calculating the Bogolubov coefficient in the near horizon limit outside the horizon for 
which T > 0. The Bogolubov coefficient is given by 

r oo 

/W = ~— / dR^R^l^^J) (113) 
Gdh JR eh 

Here oujjGdh replaces the 2u> in the standard definition of the Bogoliubov coefficient (see below for further 
explanation). Subscripts og/ig stand for outgoing/ingoing modes and R e h is the area radius at the horizon. 
Now Jgim is the metric coefficient on the space of metrics. From the form of the Hamiltonian constraint 
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this metric is seen to be diagonal in the (r, R) coordinates. In the (T, R) coordinates its form can be 
determined by using 



dr\ 2 UR x2 



9m (T,R) =gTT ^-) +SRR [ 7m ) (114) 
which yields y/gRR = 1/ T . For an expanding dust cloud the ingoing wave function is given by 

*+[r,i?,r]=eii'* r W(^+ T ). (115) 
Similarly the outgoing mode is given by 

9-[T,R,r] = et/^Wd-m/^)^ (n6) 
On the lattice the above modes become 

< = i™n e ^ (¥+TOa,i ( ii7 ) 

i 



^og = J™H e " ( 118 ) 



Here we have replaced cq by b and we have defined noj = F' , which implies T = tu/G d . With everything set 
we now evaluate (|113p by substituting for the modes (in the following we will drop the subscript i always 
remembering that we are working on a lattice). 

^=™ev(*[hu/- u ) + T(u/ + U )]) I™ d^Iexpf-^ / dR^EZ) . (119) 

G d h \G d h 2 J jR eh J- V G d h J J- J 

As a first step we calculate the Bogoliubov coefficient in the absence of the cosmological constant in which 
case J- = 1 — F/R n ~ l . We do a near horizon calculation noting that at the horizon 1 — F/R n = 0. To 
simplify the integration we define 

^y- ~ I- (120) 
With this the integral term in the above equation becomes 

1 n+l / 1 2 \ 

2F^ f°° (l + s)^r / . 4auF^ f s (l + s)^ 



, ds — k exp — i— / ds — k . (121) 

n - 1 Jo s 2 + 2s y \ G d h(n - 1) J s 2 + 2s ) y ' 

To regularize the integral at infinity we multiply it with e~ ps , p > and take the limit p — > after doing 
the integration. Thus there is a suppression for large values of s and we can assume that s is small. Thus 
retaining only the first order terms in s the integral becomes 



i 



lim / ds[s G d n n -i) e G d ft ( n-i)2 e - P s n 22) 

p^o n - 1 ' 1 ' 



The above integral evaluates to 



n — 1 



iaujF"- 1 s 
H(n - l) 2 



e G d a ( „-i) r _ (123) 

I G d /i(n - 1) I 
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With this the Bogoliubov coefficient turns out to be 

»2 ? ^F 1 / n - 1 
G d h(n-1) 



/W = ~— -exp I ~— -(w -w) +T(£j +w) I 

G d hn-1 \G d h2 J 



i2aojF 1 ' n - 1 



G d h{n - 1)2 



e o d »(n-i) r I i . (124) 

^ G d ft(n - 1) y 

Using |T(iy)| 2 = 7r/y sinh(7ry), (y real), and noting that particle creation rate corresponds to the absolute 
square of the Bogoliubov coefficient we find that 

/W = 5TT • ( 125 ) 



e G d R(n-l) _ ]_ 



By defining (To; = G d AE, AE being the energy of a shell, this can be interpreted as an eternal Schwarzschild 
black hole in equilibrium with a thermal bath at the Hawking temperature 

k B T H = ^^ T . (126) 
For instance in 3 + 1 dimensions where n = 2 we find 

k " T « = sJm' < 127 » 

We now calculate the Bogoliubov coefficient in the presence of the negative cosmological constant. We 
again work in the near horizon approximation. We note that the horizon in (jl 19)) is given by the root of T 

F 2AR 2 , . 

H n 1 n{n + 1) 

Since this equation cannot be solved generally for arbitrary n, we write the expression for T in terms of the 
surface gravity. The surface gravity is given by (for an exact expression of surface gravity see the following 
section) 

K = -2^ R = R - (129) 
However, goo = T which, as seen above, equals zero at the horizon and therefore we can approximate T 
near the horizon by 

T = —2k(R - R e h)- (130) 

Defining R — R e h := s we have that J- = —2ks and with this the expression for Bogoliubov coefficient 
becomes 

G d h \G d h 2 J Jo \2ksJ V G d h J J 

To regularize the integral we multiply it with s q e~ ps , (p,q) > (and take the limit (p,q) — ► after doing 
the integral) which would imply that we will get contributions only for those values of s which are near 
zero. For this reason we retain only the first order terms in the integral in the exponential. With this the 
above integral becomes 



= lim 



nij] f in b \ f°° / 1 i jg" n \ 

o ^-exp ( [ (u>'-u) + T(u/ + u)] ) ds(s + ^e^s^ s ) 
2nG d n \G d n * J Jo \ ) 



(132) 
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This finally evaluates to 



" r (-S' (133) 

As before we are interested in the particle creation rate for which we take the absolute square of the above 
equation which gives 

ia^GstVM— )■ (134) 

\2G d hK/ \ e G d h K _lJ 

As before, using the correspondence 010 = GdAE, we have a black hole in equilibrium with a thermal bath 
at Hawking temperature 

k B T H = — . (135) 

Since it is difficult to solve for the horizon radius for an arbitrary n we cannot give a general formula for 
the surface gravity and hence for the Hawking temperature which will show its dependence on the number 
of dimensions, the mass of the black hole and the value of the cosmological constant. However, in the 
next section we see how these three quantities control the behaviour of the Hawking temperature and the 
specific heat. 

7 Surface Gravity 

As follows from ([6]) the event horizon is defined by the condition 

F s 2Ax 2 

1 + — = 0. (136) 

x n 1 n[n + 1) 

This implies that the radius of the event horizon is given as the solution to the following equation 

2Ax n+1 + n(n + l)^™" 1 - n(n + l)F a = 0. (137) 

In particular, for n = 2 we find that the event horizon is given by 

1 (3GM A 2 + VA 3 + 9G 2 M 2 A 4 ) 1/3 /100 , 

x e h = ; — r H : • (138) 

(3GMA 2 + VA 3 + 9G 2 M 2 A 4 )V3 A v ; 

Area of the event horizon is nothing but the volume of an n-sphere of radius x e h which is given by 

n+l 

Aeh = ^rfek (139) 

where 2it~^~ /T(^-^-) is the volume of a unit n-sphere. The surface gravity is given by 

ld 9oo, ri/m s 
« = -2-^-1— ( 14 °) 

and works out to 

1 l(„ _ 4At.l 1 

(141) 



1 

K= 2 



(n - l)F a 4Ax eh 



X eh n ( n + 1)J ' 



Using F s = (1 + 2Ax 2 eh /n(n + l))x n e ^ we get 

1 



K= 2 



n-1 2Ax eh 



(142) 
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Noting that the Hawking temperature of the quantized black hole is determined by surface gravity, we 
see that this expression for surface gravity has remarkable implications for black hole thermodynamics. 
The Hawking temperature is controlled by three parameters: mass of the black hole, number of spatial 
dimensions and the cosmological constant. The presence of A could make the specific heat of the black 
hole positive (as it does for the BTZ black hole). To this effect we calculate the derivative dn/dM for fixed 
n and A. From (|142p we get that 

n-1 2A~ 



^ dn _ dx 
dM ~ dM 



2 + — 

x A n 



(143) 



where to simplify notation we write x instead of x e h. As is expected, and also evident from (|137p . x 
increases monotonically with M. Thus the condition for positivity of the specific heat is 

2Ax 2 > n(n - 1) (144) 

which is a constraint on the dimensionless quantity constructed from the cosmological constant and the 
radius of the event horizon. The above relation was also obtained in [18]. After expressing x in terms of 
M, A and n, this becomes a condition on the mutual relation between the three free parameters. A few 
things can be read off immediately: (i) if A is zero, specific heat is necessarily negative, (ii) for the 4-d AdS 
case, i.e. n = 2, specific heat is positive if Ax 2 > 1. Using the expression for x from (|138h this appears to 
translate into a complicated relation between M and A, (iii) given any n > 1, and some value of A, one 
can always choose an M, and hence an x, sufficiently large, so that the specific heat becomes positive. 



8 Conclusions 

We have used the canonical formalism to study the quantization of a spherically symmetric dust cloud in 
the presence of a negative cosmological constant in arbitrary number of spatial dimensions. It is remarkable 
that the canonical form is same in all dimensions, and is also not affected by the inclusion of the cosmological 
constant. We obtained the Wheeler-DeWitt equation which was solved to obtain the exact quantum states 
using lattice regularization. In the limit of small dust perturbations over the AdS-Schwarzschild geometry 
these states correspond to Hawking radiation. This provides justification for the choice of the inner product 
on the space of metrics. We also saw that there is an interesting interplay between the number of spatial 
dimensions, the cosmological constant and the mass of the black hole, which can render the specific heat 
of the black hole positive. 

It would also be interesting to see the role of this interplay in determining the entropy of an AdS- 
Schwarzschild black hole in arbitrary number of dimensions. This aspect will be discussed in the next 
paper where we calculate the entropy by counting the microscopic degrees of freedom using the quantized 
states obtained here. The change of the sign of the specific heat is likely to have important consequences 
for the statistical description of the thermodynamics. 
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